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BaCo2V208 is a nice example of a quasi-one-dimensional quantum spin system that can be de- 
scribed in terms of Tomonaga-Luttinger liquid physics. This is explored in the present study where 
the magnetic field-temperature phase diagram is thoroughly established up to 12 T using single- 
crystal neutron diffraction. The transition from the Neel phase to the incommensurate longitudinal 
spin density wave (LSDW) phase through a first-order transition, as well as the critical exponents 
associated with the paramagnetic-Neel phase transition, and the magnetic order both in the Neel 
and in the LSDW phase are determined, thus providing a stringent test for the theory. 

PACS numbers: 75.10.Pq, 75.50.Ee, 75.25.-j, 75.40.Cx, 75.30.Kz 



I. INTRODUCTION 

One dimensional (ID) antiferromagnctic spin chains 
have attracted a lot of attention in the last decades, 
both theoretically and experimentally. A strong quantum 
character of these systems implies that they can show 
many fascinating phenomena, not existing in systems of 
higher dimension. Also, many magnetic properties can 
be exactly calculated, owing to their low dimensionality. 
More recently, magnetic field effects have been particu- 
larly focused on, since the application of a magnetic field 
can bring new behaviors, such as magnetization plateaus, 
Bose-Einstein condensation, or novel field-induced mag- 
netic orderings. 1 Such novel phenomena were indeed 
observed respectively in e.g. the following systems of 
weakly-coupled spin chains: azurite containing distorted 
diamond spin-1/2 chains, 2 ' 3 Cu(N03)2-2.5H 2 contain- 
ing alternating spin-1/2 chains, 1 and NDMAP containing 
anisotropic Haldanc spin-1 chains. ' A great deal of exper- 
imental work has been performed on spin-1/2 and spin-1 
quantum ID systems, but very little on systems of higher 
spin. 

The quasi-lD spin-3/2 screw chain antiferromagnet 
BaCo2V2 0g is particularly interesting since it possesses a 
strong anisotropy, some sizable frustration, and, in spite 
of a larger spin, strong quantum effects. Namely, the 
low-temperature magnetic state of the Co 2+ ion (3d 7 ) in 



a distorted octahedral environment, as it is the case in 
this compound, is described by a highly anisotropic ef- 
fective spin S = 1/2. ! As a result, quantum fluctuations 
play an important role in this compound. 



In zero magnetic field, BaCo2V20s exhibits upon cool- 
ing a transition from a ID magnetic behavior to a 3D an- 
tiferromagnctic ordering at T ~ 5 K, as first evidenced 
by He et al. ' from magnetic and specific heat measure- 
ments. This system also presents a strong magnetic 
anisotropy, as shown from a twice larger susceptibility 
when the field is applied along the c— axis (chain direc- 
tion), as compared to the one in a field applied perpen- 
dicular to the c— axis. 7,8 All these results indicate that 
the chain c— axis is the magnetic easy axis in this com- 
pound. From magnetization measurements performed in 
a field applied along the c— axis, a field-induced order-to- 
disorder transition occurring above 1.8 K (critical field 
H c ~ 3.8 T at 1.8 K), instead of a spin- flop transition, 
was then evidenced,' 1 before the saturation is reached at 
H s = 22.7 T." 1 Kimura et al. 10 ' 11 attributed this order- 
to-disorder transition to the quantum phase transition 
characteristic of the S = 1/2 ID XXZ antiferromagnet 
with Ising-like anisotropy model, whose Hamiltonian is 
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given by: 

i i 

(i) 

where Sf (with a = x, y or z) denotes the a-th compo- 
nent of the spin residing at the site i of the chain, J > 
is the antiferromagnetic exchange interaction along the 
chain, fiB is the Bohr magneton, H is the longitudinal 
applied magnetic field (H || c— axis), g is the electron g- 
factor along the c— axis, and e is an anisotropic parameter 
(e < 1 for the Ising-like case). The authors obtained the 
following parameters from the fit of the magnetization 
curves measured up to 55 T for H \\ c: J/ks = 65 K (ks 
is the Boltzmann constant), e ~ 0.46, and g = 6. 2. 10,11 

An important concept of ID quantum spin systems 
in their gapless phase is the Tomonaga-Luttinger liquid 
(TLL). 1 - The transverse and longitudinal spin-spin cor- 
relation functions decay with the distance r between the 
spins asymptotically as power laws, so that they depend 
on the TLL exponent rj as follows: 

(^-(-lrirr" (2) 

and 

{S£S*} - (S z ) 2 ~ cos(2fc F r) M" 1 / 1 ', (3) 

where (S z ) is the z— component of the spin due to the 
magnetic field H applied parallel to the c— axis (0 < 
(S z ) < 1/2) and fcp is the Fermi wave number defined by 

k F = tt(1/2 - (S z )). (4) 

Mapping the spin chain onto a ID system of spinlcss 
fermions using the Jordan- Wigner transformation, 1 ' 5 this 
equation simply relates the filling of the fermion band 
1/2 — (S z ) with the wavevector 2kp connecting both 
points of the Fermi surface. Alternatively, the transverse 
component of the correlation function can be said to de- 
cay with the exponent rj x = r/ y = 77, while the longitudi- 
nal component decays with the exponent rj z = I/77. Rela- 
tion (4) implies that the value of Rf can be tuned by the 
applied magnetic field H, through {S z ). The longitudi- 
nal correlation function is thus incommensurate with the 
chain, while the transverse correlation function always 
appears at the wave number n, so that it is staggered. 

When spin chains are weakly coupled, a long-range or- 
der of the spins sets in at low enough temperature. For 
the isotropic Hcisenberg case (e = 1), the relation 77 < 1 
(i.e., t] x < Tj z ) is always satisfied in the presence of a mag- 
netic field, and thus the transverse fluctuation of Eq. (2) 
is dominant, implying a long-range ordering of the xy 
spin components (Neel order). However, for the Ising- 
like case, the relation 77 > 1 (r) x > rj z ) holds above H c 
in the low field region of the gapless phase indicating a 
dominance of the longitudinal incommensurate correla- 
tion [see Eq. (3)], before 77 < 1 (rj x < rj z ) becomes sat- 
isfied at higher field: an inversion of the TLL exponents 
rj x and rj z occurs at some critical field H* . Therefore, 



a long-range ordered phase with a longitudinal incom- 
mensurate structure is predicted theoretically between 
H c and H* before the establishment of a transverse stag- 
gered ordered phase at H > H* . Based on Eq. (3), the 
magnetic structure propagation vector of the longitudinal 
incommensurate order is expected to be: 

k = 2k F . (5) 

As {S z } varies with the field, the propagation vector fol- 
lows this variation. 

After the magnetization measurements of Hc et al. at 

1.8 K, Kimura et al. 14 performed high-field specific heat 
measurements at 200 mK, which evidenced a new mag- 
netic field-induced ordered state above H c ~ 3.9 T. The 
authors suggested that this field-induced magnetic phase 
exhibits the novel type of longitudinal incommensurate 
order described above, which is caused by the quantum 
effect inherent in the S — 1/2 quasi-lD Ising-like anti- 
fcrromagnet. This was then further explained by Oku- 
nishi and Suzuki, 1 '' who calculated the TLL exponent 77 
by the Bethe ansatz integral equations, using the values 
of J, e and g determined from the magnetization curves 
and the interchain coupling J'/ks = 0.09 K estimated 
from the mean-field theory. They found that 77 > 1 be- 
tween H c ~ 3.9 T and H* - 15 T in BaCo 2 V 2 8 .' 
This theoretical work strongly suggests that the field- 
induced transition observed below 1.8 K by Kimura et 
a/. 14 corresponds to a Neel- longitudinal incommensurate 
order transition, or, more precisely, to a transition to- 
wards an incommensurate longitudinal spin density wave 
(LSDW). Okunishi and Suzuki also calculated the transi- 
tion temperatures for the incommensurate and transverse 
staggered orders, as a function of the applied field, and 
obtained a very good agreement with the experimental 
observation for the low-field phase boundary of the lon- 
gitudinal incommensurate ordering. 1 "' 

Soon after, this novel type of field-induced incommen- 
surate order was further hinted by neutron diffraction. 16 
The magnetic 4 3 reflection was followed up to a mag- 
netic field of 5 T, and was shown to split above H c ~ 

3.9 T into incommensurate satellites at 4 3±£ . The field 
dependence of £ above H c was compared to the theoret- 
ical one but only in a narrow field range of 1 T. 

Very recently, Kawasaki et al. 1 ' determined the mag- 
netic structure of BaCo2V2 0s in zero magnetic field from 
a powder neutron diffraction experiment: the 2.18/is 
magnetic moments are aligned along the c— axis, with 
an antiferromagnetic arrangement along c and along one 
of the crystallographic axes of the square basal plane (a 
or 6), and with a ferromagnetic arrangement along the 
other axis (b or a). 

In this paper, we present a complete exploration of 
the magnetic field-temperature H-T phase diagram, up 
to 12 T and down to 50 mK, by means of single-crystal 
neutron diffraction, yielding various new and important 
results. First, besides the overall agreement of our phase 
diagram with that of Kimura et al. [i (see Fig. 11), the 
hysteretic behavior of the Neel-LSDW phase transition 
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FIG. 1: (Color online) Crystal structure of BaCo2V208. (a) Projection in the (6, c) plane: for sake of clarity, only the two 
chains with < x < 0.5 and —0.25 < y < 0.75 are represented for —0.62 < z < 1.62. (b) Projection in the (a, b) plane of the 
four chains of the unit cell (0 < x < 1, —0.25 < y < 0.75, —0.15 < z < 1.15): the arrows indicate the sense of rotation of the 
screw chains on increasing z. 



was carefully studied and a phase boundary marking 
the collapse of the LSDW phase, not predicted by the 
theoretical model of Okunishi and Suzuki, 15 was evi- 
denced at around 9 T. Second, a precise study of the 
critical exponents was performed in the low field part 
(0 < H < H c ~ 3.9 T) of the phase diagram. Third, the 
antiferromagnetic structure was determined in zero field, 
in good agreement with the one recently determined from 
powder neutron diffraction, 1 ' and the structure of the 
field-induced LSDW phase was determined for the first 
time. Last, a very detailed survey of the propagation 
vector and of the intensity of chosen magnetic reflections 
was performed up to 9 T, yielding an excellent agreement 
with the predicted field dependence of the incommensu- 
rate modulation over a large field range. 



The paper is organized as follows: after the introduc- 
tion, the crystal structure and the experimental proce- 
dures are described in Sec. II. In Sec. Ill, the results 
are then presented and analyzed. In subsection III A, a 
characterization of the sample by means of specific heat 
measurements is presented, followed by the presentation 
of the neutron scattering measurements. The determi- 
nation of the nuclear and magnetic structures at H = 
are presented in subsections IIIB and IIIC, respectively 
In subsection HID, the temperature and magnetic field 
scans are presented, including a study of the critical ex- 
ponents and a careful survey of the incommensurate mag- 
netic peak position as a function of the magnetic field, 
yielding the determination of the H-T phase diagram. 
In subsection HIE, the magnetic structure in the LSDW 
phase is presented. Sec. IV is devoted to a discussion, 
followed by a conclusion in Sec. V. 



II. EXPERIMENTAL 
A. Crystal structure and growth 

BaCo2V2C>8 crystallizes in the centrosymmetric tetrag- 
onal body-centered IAi/acd (No. 142) space group, with 
a = 12.444 A, c = 8.415 A, and 8 chemical formulae 
per unit cell. 18 The sixteen magnetic Co 2+ atoms of the 
unit cell are equivalent (Wyckoff site 16/). The spin- 
3/2 Co 2+ atoms are arranged in edge-sharing CoOg oc- 
tahedra forming screw chains, running along the c— axis, 
and separated by non-magnetic V 5+ and Ba 2+ ions (see 
Fig. 1). As schematized by the arrows in Fig. 1(b), a 
screw axis 4i runs through the first type of chains of 
the lattice, transforming Coi (Coy) at z = 1/8 into C02 
(C02') at z = 3/8, C03 (C03') at z = 5/8, and then C04 
(C04/) at z = 7/8 (anti-clockwise rotation as z increases), 
while a screw axis 43 runs through the second type of 
chain of the lattice, transforming C05 (C05/) at z = 1/8 
into C06 (C06') at z = 3/8, C07 (C07/) at z = 5/8, and 
then Cos (Cog/) at z = 7/8 (clockwise rotation as z in- 
creases). The C01/...C04' and Co5<...Cos' chains are ob- 
tained from the C01...C04 and C05...C08 chains, respec- 
tively, simply by applying the /—centering of the lattice: 
Coi at z — 1/8 yields C03' at z — 5/8, Co 5 at z = 1/8 
yields C07' at z = 5/8, by translation of (-5, 5, -|). 

Two BaCo2V2 08 single-crystals were grown at Institut 
Neel (Grenoble, France) by the floating zone method. 19 
In both cases, a 5 cm long cylindrical crystal rod, of 
3 mm diameter, was obtained, with the growth axis at 
about 60° from the c— axis. From the first rod, a 2.2 mm 
thick slice was cut perpendicular to the c— axis, yielding a 
135 mm 3 sample with elliptical faces of about 3 x 6.5 mm 2 
(first crystal devoted to the neutron diffraction experi- 
ments). Nearby on the same rod, a very small crystal of 



4 



2.5 mg was cut with faces perpendicular to the c— axis for 
the specific heat measurements. From the second rod, the 
sample was cut in order to have a cubic shape, with about 
3 mm long edges along the crystallographic axes (second 
crystal also devoted to the neutron diffraction experi- 
ments). Note that this second neutron sample, though 
5 times smaller, was much better suited for structure re- 
finements, in particular because of the strong neutron 
absorption by Co: the various neutron paths (and thus 
the absorption) were then very similar for all measured 
reflections. However, the first sample with a much larger 
volume was more appropriate for the measurements of 
the incommensurate satellite peaks with small intensity. 



B. Specific heat measurements 

The specific heat measurements were performed at 
CEA-Grenoble / SPSMS-IMAPEC (Grenoble, France) 
on a commercial Physical Properties Measurement Sys- 
tem (PPMS, Quantum Design) equipped with a 3 Hc in- 
sert using a semi-adiabatic heat-pulse technique. On this 
apparatus, magnetic fields up to 9 T could be reached and 
the temperature was varied between ~ 600 mK and 7 K. 
The 2.5 mg sample was mounted with its cutting surface 
parallel to the sample platform and thus with the c— axis 
vertical, parallel to the applied magnetic field. 



C. Neutron diffraction experimental procedures 

All neutron diffraction experiments were performed at 
the Institut Laue Langevin (ILL) high-flux reactor in 
Grenoble, France. The two single- crystals were previ- 
ously aligned with the c— axis vertical on the OrientEx- 
press neutron Laue diffractometer. A series of three ex- 
periments were performed on the first single-crystal (with 
an elliptical shape). The first one was performed on the 
CSIC/CEA-CRG D15 single-crystal diffractometer, op- 
erated in the 4-circle mode, with the sample mounted in 
a 2 K - displex refrigerator. Two experiments were then 
performed on the CEA-CRG D23 single-crystal two-axis 
diffractometer with a lifting arm detector. The sample 
was mounted in a 6 T then a 12 T vertical field cryo- 
magnet, equipped with a dilution insert, with the field 
applied along the c— axis. Later on, a fourth experi- 
ment was performed on D23 on the second (cubic shaped) 
single-crystal (not available at the time of the previous 
experiments). The sample was placed in a standard ILL 
orange cryostat with the c— axis vertical. Incident wave- 
lengths of 1.173 A and 1.280 A were used, respectively 
on D15 and D23. 

For the zero-field experiments on D15 and D23, no 
collimation was used: only diaphragms were put before 
and after the sample, in order to optimize the signal over 
noise ratio. As concerns the experiments performed on 
D23 with the 6 and 12 T cryomagnets, a 20' vertical 
collimation was put before the sample, in addition to 



diaphragms, in order to improve the vertical resolution 
(along the reciprocal lattice vector c*) and thus to sep- 
arate the satellites appearing in the LSDW phase at ±£ 
(£ << 1 close to H c ). Doing so, a diaphragm of verti- 
cal aperture 4 mm could be installed before the detector, 
without cutting any signal, yielding, together with the 
use of the collimation, a huge reduction of background. 
A resolution of about 0.07 r.l.u. (reciprocal lattice units) 
could then be achieved along the c* direction. Putting 
the c— axis of the sample vertical on D23 implied that hkl 
reflections up to I — ±3 could be collected with the or- 
ange cryostat, while Q = (Qh, Qk, Ql) scattering vec- 
tors with -0.5 < Q L < 2.15 and -0.3 < Q L < 1.1 could 
be reached with the 6 T and 12 T cryomagnets, respec- 
tively. 

The zero-field experiment on D15 performed on the 
first single-crystal was aimed at determining the mag- 
netic structure in the Neel phase at T = 2 K, together 
with the nuclear one above and below the transition, and 
at doing a precise temperature dependence of magnetic 
peaks. This structure determination was later improved 
on D23 using the second single-crystal mounted in a cryo- 
stat. The experiment done on D23 with the 6 T cryomag- 
net was aimed at refining the nuclear and magnetic struc- 
tures in the LSDW phase at H = 4.2 T and T = 50 mK 
and at measuring temperature and field dependences up 
to 6 T. These field and temperature scans were then com- 
pleted at higher field using the 12 T cryomagnet on D23. 



III. RESULTS AND ANALYSIS 
A. Specific heat measurements 

Before performing any neutron diffraction measure- 
ment, our BaCo2V2 08 single-crystal was first character- 
ized by X-rays 10 and specific heat. The aim of the specific 
heat measurements, performed in a magnetic field ap- 
plied along the c— axis, was to obtain a macroscopic sig- 
nature of the magnetic-paramagnetic phase boundaries 
under a magnetic field. The specific heat was measured 
as a function of the temperature for various magnetic 
field values ranging between and 9 T (see Fig. 2). 

In zero magnetic field, a sharp A— type anomaly is seen 
at T c ~ 5.5 K evidencing a second-order phase transition 
between the paramagnetic and Neel phases when lower- 
ing the temperature. The critical temperature T c and 
the amplitude of the A— like peak both decrease across 
the paramagnetic-Neel phase transition as H increases. 
A slower decrease of T c is observed for H > 4 T, i.e., for 
the paramagnetic-LSDW phase transition. For H > 7 T, 
the transition is no longer visible (the measurements per- 
formed at 8 and 9 T are not shown), either because the 
critical temperature is lower than 600 mK (minimal tem- 
perature reached during the measurement) or because the 
anomaly is too weak to be observed. From all the mea- 
surements shown in Fig. 2, the ordcr-to-disordcr transi- 
tion line in the H-T phase diagram was obtained (see 
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FIG. 2: (Color online) Temperature dependence of the specific 
heat measured in BaCo2 V2O8 in different magnetic fields from 
to 7 T applied along the c— axis. 

Fig. 11) prior to the neutron diffraction measurements. 



B. Nuclear structure at H — 

The first step of the neutron diffraction experiments 
consisted in refining the low-temperature nuclear struc- 
ture, previously determined from X-ray diffraction at 
room temperature. ltS This was done both in the para- 
magnetic phase and in the Neel phase, in order to check 
that no modification of the nuclear structure occurs when 
entering the Neel phase. To do so, more than 300 inde- 
pendent reflections, allowed in the Mi/acd space group, 
were collected both at T = 14 K and T = 2 K on D15 
using the first crystal. About a hundred independent for- 
bidden reflections were also measured and checked to be 
null. Though the weighted least-squares R— factor was 
not completely satisfactory (Rf^w ~ 10%), due to the 
fact that cobalt absorbs neutrons a lot and that this sam- 
ple has a very asymmetric shape (dimensions of about 
6.5 x 3 x 2.2 mm 3 ), these refinements allowed to conclude 
that the nuclear structure remains identical across the 
magnetic phase transition, with no change in the atomic 
coordinates, within the error bars. 

We now present in detail the experiment performed on 
D23 at T = 1.8 K on the second crystal, which yields a 
much better refinement owing to its cubic shape. The fol- 
lowing lattice parameters were obtained, from centering 
22 reflections: a = 12.404(2) A and c = 8.375(12) A. 622 
nuclear reflections hkl up to I = ±3 (among them 172 
independent ones), allowed in the IAi/acd space group, 
were then collected. About 50 independent forbidden re- 
flections hkl with h + k + I even (i.e., nuclear reflections 
allowed by the /—centering of the lattice but forbidden 
either by the screw axis 4i or by a glide plane of the space 
group) were measured as well and checked to be null. A 
few hkl reflections with h + k + l odd were also measured 
at large Q (in order to have a null magnetic contribu- 



TABLE I: Atomic coordinates and isotropic Debye-Waller 
factors obtained from the nuclear structure refinement per- 
formed on D23 on the cubic shaped crystal at T = 1.8 K 
in I4i/acd space group (agreement Rp — factor = 2.38% and 
R F 2 W — factor = 4.75%). The larger error bars for the z co- 
ordinates come from the fact that the c— axis was set vertical 
on the diffractometer, limiting the I Miller index to ±3. 
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tion) and their zero intensity confirmed the /—centering 
of the lattice. The nuclear refinement was done with the 
FULLPROF software. 2 " In addition to the scale factor, the 
x,y,z coordinates and the Bi so isotropic Debye- Waller 
factors were refined for all atoms except vanadium, 21 to- 
gether with the A/2 ratio (smaller than 0.1%) and the ex- 
tinction parameters. 22 As expected, a much better agree- 
ment Rp2 w — factor than for the first crystal was obtained 
(Rf 2 w = 4.75%), as a consequence of a better suited crys- 
tal shape, yielding the atomic coordinates and isotropic 
Debye-Waller factors listed in Table I. Let us emphasize 
that the atomic coordinates are very similar to those de- 
termined on the first crystal and that the Debye-Waller 
factors have reasonable values. Last, the four refined ex- 
tinction parameters are smaller than the six ones refined 
on the first crystal, which is consistent with the smaller 
crystal size. The absorption corrections of the intensities 
(using the AVEXAR software from the CCSL library 23 ) 
did not improve the refinement and was not done in the 
following magnetic refinements. 

C. Magnetic structure in the Neel phase (H = 0) 

The first observation by Kimura et al. of the mag- 
netic signal on the 4 3 reflection allows to easily deduce 
the propagation vector of the magnetic structure in the 
Neel phase. The magnetic signal appears at H = at 
positions hkl with h + k + l odd integer, that is, on the 
nuclear peaks forbidden by the body-centering / of the 
lattice. The propagation vector of the magnetic structure 
is thus \<laf = (1, 0, 0) [or, equivalcntly, k^F = (0, 1, 0)] 
owing to the fact that a > c. > The magnetic structure 
was first determined on D15 using the crystal of ellipti- 
cal shape and then on D23 using the cubic one. Only 
the latter experiment will be presented again, since both 
gave the same result with a better agreement factor for 
the cubic shaped crystal. 313 different magnetic reflec- 
tions, reducing to 143 independent ones, were collected 
on D23 on the second crystal, counting 10 seconds per 
point. The scale factor, the x atomic coordinate and 
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the isotropic Dcbye- Waller factor of cobalt, the four ex- 
tinction parameters, and the A/2 ratio obtained from the 
nuclear structure refinement provided us with all the nec- 
essary parameters for a quantitative analysis of the mag- 
netic structure. 

For a second-order phase transition (as it is the case 
here) from the paramagnetic phase to the Neel phase, the 
group theory can be used in order to predict the possible 
spin structures compatible both with the space group of 
the nuclear structure and with the propagation vector of 
the magnetic structure, according to the irreducible rep- 
resentations (irreps) t v , as described in Ref. 26. As a re- 
sult, this usually reduces the number of independent pa- 
rameters to be refined. The BASIREPS software from the 
FlJLLPROF suite was used and four irreducible represen- 
tations of dimension 2 were found, occurring 3 times each 
in the group representation: T = 3t 1 ©3t 2 ©3t 3 ©3t 4 . In 
the present case, the use of group theory thus reduces the 
number of fitting parameters from 24 down to 6, what- 
ever the irreducible representation. The four possible 
magnetic arrangements, deduced from this analysis, are 
summarized in Table II. Note that the Coi'...Cos' atoms 
do not appear in this table. The magnetic moments of 
the primed cobalt atoms are simply obtained from the 
unprimed ones as follows: the cobalt atoms that are re- 
lated by the / lattice translation (5,5,5) have antiparal- 
lel magnetic moments, owing to the propagation vector. 

TABLE II: (Color online) Magnetic components for the eight 
Co atoms: Co! (0.331, -0.081, 1/8), Co 2 (0.331, 0.081, 3/8), 
Co 3 (0.169, 0.081, 5/8), C04 (0.169, -0.081, 7/8), Co 5 (0.169, 
0.419, 1/8), Co 6 (0.169, 0.581, 3/8), Co 7 (0.331, 0.581, 5/8), 
and Cos (0.331, 0.419, 7/8), for each irreducible representa- 
tion t , t 2 , t 3 , and t 4 . The non-circled signs ('+' and ' — ') 
and the circled signs ('©' and 'B') correspond to independent 
magnetic moments. 
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Looking at Table II shows that either the two types of 
chains (C01...C04 and C05...C08) or two successive planes 
along c (at a distance of Az = 1/4) are magnetically in- 
dependent, depending on the irreducible representation, 
r 1 and r 2 or t 3 and r 4 , respectively. Six parameters 
are thus to be refined in all cases: the magnetic compo- 



nents (m x , m y , m z ) of C01 and {m' x , m' y , m' z ) of C05 or 
C02, depending on the irreducible representations we are 
dealing with, r 1 and r 2 or r 3 and r 4 , respectively. 

The refinement of the magnetic structure, detailed in 
the following, was tested for each irreducible represen- 
tation, and only r 1 gave a satisfactory agreement, with 
zero m x ,m' x ,my, and m' y (within the error bars). This 
confirms the prediction of Ising-likc chains, with the easy 
axis along c, and that the exchange couplings are an- 
tiferromagnetic along the screw chains. The magnetic 
components along x and y were then fixed to zero. 




50 100 150 200 

Observed Intensity (arb. units) 



FIG. 3: (Color online) Graphical representation of the mag- 
netic structure refinement in the Neel phase of BaCo2V2 0s at 
T = 1.8K and H — 0: calculated versus measured integrated 
intensities. 

Note that Kimura et al. 16 claimed that the absence of 
magnetic signal on the 01 position evidences that the 
magnetic moments are aligned along the c— axis, due to 
the fact that neutrons are sensitive only to the trans- 
verse component (that is, perpendicular to the scatter- 
ing vector) of the magnetic moment. But this absence 
of intensity also occurs for peaks 100 and 010, so that 
no conclusion can be given on the magnetic moment di- 
rection from this single observation: the zero intensity in 
this case rather comes from the magnetic structure factor 
(magnetic arrangement within the unit cell) . A complete 
refinement was thus necessary, as done in the present 
work, in order to conclude on the moment orientation. 

The stabilized magnetic structure with the Ising char- 
acter has a lower symmetry than the paramagnetic space 
group, and in particular is no longer invariant under the 
symmetry operator — y + j,x — j, z + j (4i screw axis). 
Two magnetic domains are thus expected, related by this 
symmetry element, with m z = m' z for the first domain 
and m z = —m' z for the second domain, imposing the 
same value of the magnetic moment on both types of 
chains. The refinement of the magnetic structure was 
then performed, for the r 1 irreducible representation, re- 
fining m z and the domain population. The following re- 



TABLE III: Observed intensities 7 b s of a few magnetic reflec- 
tions hkl, ordered by increasing sin#/A, compared to the cal- 
culated ones: Z^aic and /jfaic are the contributions of each mag- 
netic domain (with respective populations 49.5% and 50.5%) 
to the total intensity J*"^ = J^ alc + Icaic- Note that for some 
of them, only one domain contributes, e.g., for 2 3 and 3 2 0, 
while for some others, both domains contribute. The small 
difference in the total intensity between two reflections hkl 
and khl is due to the small difference between the domain 
populations. 
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suits were found, with a weighted least-squares R— factor 
of 11.3% (R F = 8.3%): m z = 2.267(3)^ s /Co 2+ with 
domain populations of 49.5(4)% and 50.5(4)%, that is, 
equally populated domains, as expected. This moment 
amplitude of about 2.3/is/Co 2+ is consistent with a mag- 
netic moment of reduced amplitude, as compared to the 
expected one for spin-3/2 (m = 3/jb), in presence of 
quantum effects. The calculated intensities plotted as 
a function of the observed ones emphasize the quality 
of the fit (see Fig. 3). Note that refining the magnetic 
structure with a single domain and with no constraint 
between both types of chains yields, with the same agree- 
ment factor and the same calculated intensities, two very 
different magnetic moments on both chains (3.206 and 
0.016^b/Co 2+ ), which is not realistic. 

The calculated and observed intensities of a few mag- 
netic peaks are listed in Table III. Looking, e.g., at re- 
flections 2 3 and 3 2 shows unambiguously that the 
presence of both domains is mandatory to explain our 
data. Figure 4 shows the magnetic structure in both mag- 
netic domains. In the first domain, the C01...C04 chain 
has magnetic moments parallel to those in the C05...C08 
chain and antiparallel to those in the Co5'...Cos' chain, 
while it is the contrary in the second domain. One can no- 
tice that the 4i screw axis running through the C01...C04 
chain transforms the C01-C07' pair into the C02-C08 
pair [see Fig. 4(b)]. Identical interactions couple both 
pairs which however have a different spin arrangement 
(antiparallel for C01-C07' and parallel for C02-C08, in 
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Magnetic domain #1 




Magnetic domain #2 




FIG. 4: (Color online) Magnetic structure of BaCo2V20g in 
the two magnetic domains of the Neel phase (H = and 
T — 1.8 K): perspective view (left figures) and projection 
along the c-axis (right figures). The cobalt atoms and their 
magnetic moment are labeled as in Fig. 1. The two kinds of 
chains are plotted in different colors: an axis 4i runs through 
the red chains while an axis 43 runs through the blue ones. 
Looking at the figures on the right allows an easy visualization 
of the connection between the two magnetic domains: the 
application, for example, of the 4i screw axis running along 
the C01...C04 chain transforms domain #1 into domain #2 
(C01 transforms into C02 while Co 7 / transforms into Cos). 



the first domain). This suggests the presence of sizable 
frustration between the two kinds of magnetic chains, 
which is further argued in Ref. 27, owing to a detailed 
analysis of the exchange couplings. 



Last, our results compare well with the very recent 
ones of Kawasaki et al. 1 *' These authors refined the same 
magnetic structure as the one we obtained for domain 
#2, using neutron powder diffraction, which is not sensi- 
tive to the domains. Their refined value of 2.18/is/Co 2+ 
for the magnetic moment is slightly smaller than the 
2.267(3)/is/Co 2+ obtained in the present study. 
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D. Field and temperature dependences, critical 
exponents, and H-T phase diagram 

In this subsection, we present the magnetic field and 
temperature dependences of selected magnetic reflec- 
tions, in order to obtain the critical exponents, the field 
dependence of the propagation vector and of the mag- 
netic intensities, and the phase diagram that is summa- 
rized in Fig. 11. Let us recall that when going from 
the Neel to the longitudinal incommensurate magnetic 
phase (LSDW phase), the antiferromagnetic peak splits 
above H c into two incommensurate satellites at posi- 
tions ±£ on either side of the previous peak along c* 
(see introduction). Therefore, the propagation vector 
k-AF = (1, 0, 0) in the Neel phase becomes, by conti- 
nuity, k-LSDW — (li 0) Oi so that the incommensurate 
satellites are observed at H ± \s.lsdw with H a nuclear 
scattering vector of the I lattice. 28 

Figure 5 (a) presents some temperature scans per- 
formed across two different phase transitions: (i) a se- 
ries of temperature scans performed on top of strong 
2 3 magnetic peak, across the Neel-paramagnetic phase 
transition, for different field values between and 3.25 T 
(below H c ), (ii) a temperature dependence of the 2 3£ re- 
flection with £ = 0.127, across the LSDW-paramagnetic 
phase transition, at H = 5 T (above H c ). The latter 
curve was obtained from Ql — scans performed across one 
of the two incommensurate magnetic satellites ±£: fitting 
all this series of scans by a Gaussian function allowed to 
conclude that the position of the peak does not change 
at all with temperature and to obtain the temperature 
dependence of the height of the peak. The transition 
from the magnetic phases to the paramagnetic one can 
be clearly seen from the vanishing of the probed reflec- 
tions. Note that in the LSDW phase, the magnetic sig- 
nal is much weaker (more than 10 times smaller) than in 
the Neel phase. Last, it is worth noting that the tem- 
perature scans performed in zero field for increasing and 
decreasing temperatures give exactly the same curve, evi- 
dencing the absence of any temperature hysteresis for the 
Neel-paramagnetic phase transition and thus in favor of 
a second-order phase transition. 

These measurements allowed to determine the critical 
exponent ft for the temperature-induced phase transition 
from the long-range antiferromagnetic ordering (Neel 
phase) to the paramagnetic phase for different magnetic 
fields."" The order parameter (i.e., antiferromagnetic mo- 
ment m) can be expressed as follows: m(T) oc (T c — T)' 3 , 
so that the neutron intensity curves were fitted to the ex- 
pression I(T) = I +A(1 - T/T c ) 2fi , where I and A are 
the background intensity and the magnetic intensity at 
T = 0, respectively. In our fitting procedure, the back- 
ground intensity io was first determined from the high 
temperature data and then A, ft and T c were treated as 
the fitting parameters. The method, described hereafter, 
that was used to extract the values of the critical expo- 
nent ft and the critical temperature T c , is the same as the 
one described e.g. in Ref. 30. As exemplified in Figs. 5(b) 
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FIG. 5: (Color online) (a) Temperature dependence of the 
magnetic signal: neutron counts on top of the strong 2 3 
reflection of the Neel phase for various magnetic field values 
between and 3.25 T and on top of the 2 3 0.127 reflection 
of the LSDW phase at H = 5 T. The latter curve was ob- 
tained from Ql— scans performed across the 2 3 0.127 position 
counting one minute per point (and then renormalized to 10 
seconds for the figure). The black solid lines are power- law 
fits performed on a temperature range ST — T c — T ~ 1 K 
as explained in the text, (b) and (c) The set of the ft and T c 
values as obtained for various fitting temperature ranges ST 
at H = together with a linear extrapolation of both param- 
eters to ST — for a precise determination of their values (see 
text). Inset of panel (a): log-log plot of the magnetic inten- 
sity with subtracted background, I — Io, as a function of the 
reduced temperature, 1 — T/T c , at H — 0. The data (empty 
circles) and the power-law fit (solid line) are presented using 
the extrapolated values for ST — ^ 0: T c = 5.556 K (for the 
data and the fit) and ft = 0.307 (for the fit). 



and (c) for the case of H = 0, power-law fits to the data 
were performed over a progressively shrinking temper- 
ature range ST = T c — T, and the obtained ft and T c 
values were found to increase smoothly, roughly linearly, 
with decreasing temperature range ST. The solid lines in 
Fig. 5(a) are such fits for the various temperature scans 
performed on a temperature range ST ~ 1 K. Below some 
critical value (5T) m ; n of the fitting temperature range, ft 
and T c were found to suddenly diverge. This divergence is 
due to the rounding of the transition very close to T c that 
prevented us from reducing the fitting range too much. 
For this reason the set of ft and T c values obtained for a 
set of temperature ranges ST > (<5T) min were fitted by a 
linear regression as a function of ST, in order to extrap- 
olate the values of ft and T c for ST — > [see solid line in 
Figs. 5(b) and (c)]. These extrapolated values were con- 
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sidered as the final experimental values of /3 and T c . For 
the four different fields, no significant variation of the 
extrapolated j3 value was found (values between 0.307 
and 0.328), as expected since we are dealing with the 
same universality class. These values yield an average 
critical exponent f3 = 0.32(1), which is consistent with 
the numerical estimates for the critical exponent of the 
three-dimensional Ising universality class, 33 0.308-0.332, 
to which BaCo2V20s belongs, and with the one obtained 
in recent experimental works done in zero field. 1 '' 32 The 
critical temperature obtained in zero field by extrapola- 
tion [see Fig. 5(c)] amounts to T c = 5.556(2) K. To estab- 
lish the relation of the applied procedure for extracting 
the values of f3 and T c to the more familiar one, the zero- 
field temperature scan together with the corresponding 
power-law fit are also presented in the following way: log- 
log plot of the magnetic intensity I — Iq as a function of 
the reduced temperature 1 — T/T c [see inset of Fig. 5(a)], 
where the corresponding extrapolated values of (3 and T c 
are used. The power-law fit nicely follows the data up 
to 1 — T/T c ~ 0.1 in the log-log representation, which 
further illustrates the good quality of the fit. 




FIG. 6: (Color online) Field dependence of the magnetic sig- 
nal at various temperatures between 1.85 and 5 K: neutron 
counts on top of the strong 2 3 reflection of the Neel phase. 
The black solid lines are fit to a power law as explained in the 
text. Inset: /?' values as a function of the temperature (see 
text). The dashed black line is a linear fit of these data. 

Figure 6 presents some magnetic field scans, performed 
on top of the 2 3 reflection, for different tempera- 
tures ranging between 1.85 and 5 K, across the Neel- 
paramagnetic phase transition. Some of them were done 
both for increasing and decreasing field (not shown here), 
resulting in the absence of any field hysteresis and thus 
confirming further the prediction of a second-order Neel- 
paramagnetic phase transition (field-induced order-to- 
disorder transition of He et al. mentioned in the intro- 
duction). Here again, the same analysis as for the tem- 
perature scans was performed, now fitting the neutron 
intensity curves on a progressively shrinking field range 
5H = H c -H to: 1(H) = I + A' (1 - H/H c ) 2 f 3 ' ', with 



Iq, A' , H c , and f3' as fitting parameters. The solid lines 
in Fig. 6 correspond to such fits performed on a field 
range 5H ~ 0.6 T. The /3' values, obtained by extrapo- 
lation of the f3'(SH) curves to H — > 0, are plotted in the 
inset of Fig. 6 as a function of the temperature. This ob- 
tained exponent j3' decreases in a linear-like fashion from 
0.32(3) at T = 5 K down to 0.12(2) at T = 1.85 K and 
extrapolates to ft' ~ 0.00(2) at T = (see dashed line 
in the inset of Fig. 6). This result can be understood on 
the basis of the behavior observed in temperature scans 
described above, as will be discussed in section IV. 
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FIG. 7: (Color online) Field dependence of the magnetic sig- 
nal at T = 50 mK, measured on top of 2 30 for two different 
field sweep rates, 0.2 T/min (a) and 3 mT/min (b), with 
increasing and decreasing field. A sizable field hysteresis is 
observed for the Neel-LSDW phase transition, whose magni- 
tude increases with the sweep rate (AH ~ 0.2 T at 3 mT/min, 
~ 0.32 T at 0.2 T/min). 

Fig. 7 presents some magnetic field scans, performed 
on top of the 2 3 magnetic reflection of the Neel phase 
at T = 50 mK, across the Neel-LSDW phase transi- 
tion, for increasing and decreasing magnetic fields at two 
different sweep rates. One can notice the presence of 
a sizable field hysteresis, whose magnitude strongly de- 
pends on the field sweep rate, as well as the field de- 
pendence of the intensity between about 3 and 4 T [flat 
across the complete Neel phase for a field sweep rate of 
0.2 T/min in Fig. 7(a) but not for a slower field sweep 
rate in Fig. 7(b)] . This last result will be discussed in sec- 
tion IV. All these observations evidence the first-order 
nature of the Neel-LSDW phase transition, consistent 
with the magnetocaloric measurements of Kimura et al. 1 1 
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Such field scans were also performed at different temper- 
atures between 50 mK and 1.8 K (not shown here) in 
view of completing the phase diagram that is presented 
in Fig. 11. 
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FIG. 8: (Color online) Ql — scans performed across the Q = 
(2, 3, 0) scattering vector at T = 50 mK for different increas- 
ing magnetic field values, (a) Neel phase, (b) Coexistence of 
the Neel and LSDW phases, (c) and (d) LSDW phase. For 
a better visibility, the scans measured at 4 and 3.98 T in (b) 
were shifted vertically by 2000 and 4000, respectively. 

In order to follow the splitting of the antiferromagnetic 
peak across the Necl-LSDW phase transition, Ql — scans 
(along c*) were performed centered on the strong 230 



magnetic peak at many different field values between 
and 10 T. Figure 8 presents a few of these scans, mea- 
sured at T = 50 mK for an increasing magnetic field. 
For H < 3.92 T [see Fig. 8(a)], one single magnetic 
peak is observed in the Neel phase: this peak, at ev- 
ery field value, could be well fitted by a Gaussian func- 
tion located at Ql = with half width at half maxi- 
mum SQ L = 0.0715 r.l.u. For 3.94 < H < 4.12 T [sec 
Fig. 8(b)], two satellites are observed at Ql = ±£ in 
addition to the previous magnetic peak at Ql = 0, ev- 
idencing clearly the coexistence of the Neel and LSDW 
phases. This gives an additional strong evidence for the 
first-order nature of the Neel-LSDW phase transition. In 
this field range, all scans could be well fitted by three 
Gaussian functions, at Ql — and at Ql = ±£. These 
fits allowed to extract the positions ±£ of the satellites 
of the LSDW phase and the intensity of the three peaks, 
while their width was fixed to the SQl value determined 
from the low field scans. Note that the scans measured 
between 4.02 and 4.12 T could as well be fitted by two 
Gaussian functions, i.e., with no central peak, but with 
a larger width. For H > 4.14 T [see Fig. 8(c)], the cen- 
tral peak has completely disappeared, contrary to the 
observations of Kimura et al., while the two satellites 
progressively move away from Ql = 0. This discrep- 
ancy might come from a different quality of their crystal 
or a non negligible rate of A/2 in their measurements. 
As can be seen in Figs. 8(c) and (d), a complete scan 
was performed up to 6 T, fitted by two Gaussian func- 
tions symmetric with respect to Ql — 0, while only the 
positive side in Ql was scanned from 6.5 T (fit by one 
Gaussian function). The scan at H = 9.25 T was counted 
for 10 minutes per point and the satellite peak located at 
Ql = 0.308 r.l.u. is hardly visible. At H = 9.5 T, noth- 
ing was visible after 10 minutes of counting, indicating 
that the LSDW phase no longer exists above H p ~ 9.25 T 
at T = 50 mK. 

The previous measurements and analysis were also per- 
formed for decreasing field and allowed to obtain pre- 
cisely the magnetic field dependence of the positions Ql 
of the magnetic peaks (see Fig. 9). In addition, the field 
dependence of the integrated intensity was also deter- 
mined for the magnetic peak of the Neel phase, 2 3 0, 
and for the two satellites of the LSDW phase, 2 3 ±£ (see 
Fig. 10). 

In Fig. 9, a hysteresis of nearly 0.2 T can be seen (em- 
phasized in the inset of Fig. 9), as already observed from 
the field dependence of the signal measured on top of 
the 2 30 magnetic peak (Fig. 7). This width of nearly 
0.2 T corresponds to the field range over which both the 
Neel and the LSDW phases coexist (see inset of Fig. 9). 
The incommensurate modulation of the propagation vec- 
tor, £, continuously increases when the field increases, as 
expected; in other words, it never locks to some com- 
mensurate value. Indeed, owing to the fact that there 
are four Co 2+ ions along the chain in a unit cell and 
using Eqs. (4) and (5), it is predicted that the incom- 
mensurate modulation of the propagation vector must 
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FIG. 9: (Color online) Field dependence of the position along 
c* of the magnetic reflections (red and blue circles for increas- 
ing and decreasing field, respectively), obtained from the fit 
of the Ql— scans: antiferromagnetic peak at Ql — in the 
Neel phase, incommensurate satellites at Ql ~ ±£ in the 
LSDW phase with £ the incommensurate modulation. The 
inset presents a zoom of these curves around the transition. 
The solid black line in the main panel corresponds to 4M z c °7j, 
where M™ is the magnetization curve M Z (H) of Ref. 14 cor- 
rected for a residual Van Vleck contribution and g = 6.0 as 
described in the text. 



be proportional to the field-induced uniform magnetiza- 
tion M z along c as follows: £ = 4 (S z ) = AM Z J (gLis). 
In order to compare our results to this prediction, we ex- 
tracted the magnetization curve from Ref. 1 4 and further 
corrected it for a residual Van Vleck paramagnetic contri- 
bution of 0.00733 [ji B /T (Ref. 3) to obtain M z cor . Then 
4M™ r /p = M 2 cor /1.5 was plotted in Fig. 9 (see black solid 
line) with M™ 1 ' in fig and using g = 6.0 as determined 
from the value of M™ r at saturation (i.e., 3.0{1b/Co 2+ at 
30 T). The very good agreement between M^ OT (H) and 
Ql(H) = £(H) in Fig. 9 confirms the prediction of the 
theory for the LSDW phase in BaCo 2 V208- 

Figure 10 shows that the integrated intensity of each 
satellite (LSDW phase) is, just above the transition, 
about 5 times smaller than the one of the commensu- 
rate reflection (Neel phase) , while it decreases rapidly as 
the field increases between about 4 and 9 T, and then 
completely disappears. Both satellites have exactly the 
same intensity (within the error bars) on the whole field 
range as expected, since they are equivalent and have 
the same Q modulus and thus the same magnetic form 
factor. Here again, the hysteresis of nearly 0.2 T at the 
transition is well visible (see the inset of Fig. 10). 

Last, all the exhaustive temperature and field depen- 
dences, like those shown in Figs. 5, 6, 7, and 8, allowed us 
to obtain the critical temperatures and critical fields at 
different magnetic fields and temperatures, respectively, 
and thus to obtain a detailed H-T phase diagram pre- 
sented in Fig. 11 (red circles). For comparison, we add 
to Fig. 11 the points obtained from our specific- heat mea- 
surements as well as the points published in Ref. 14. 
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FIG. 10: (Color online) Field dependence of the integrated 
intensity of the magnetic peaks, obtained from the fit of the 
Ql— scans for increasing (closed circles) and decreasing (open 
circles) field. The inset presents a zoom of these data around 
the transition. The scattering of the red and blue data sets 
reflects the typical error bars that should be considered. 
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FIG. 11: (Color online) Magnetic H-T phase diagram of 
BaCo2V208 obtained from the specific heat measurements 
presented in subsection III A (black squares) and from the 
present neutron diffraction measurements (red circles). As 
concerns the Neel-LSDW phase transition (see zoom in the in- 
set), the closed and open red circles correspond to increasing 
and decreasing field scans, respectively. The blue open circles 
in the main panel are the data of Kimura et al. obtained 
from macroscopic measurements and the blue dashed line is 
the theoretical LSDW-paramagnetic transition line predicted 
by Okunishi and Suzuki. 15 



Except for the last two points at 9 and 9.25 T, our 
phase diagram is in very good agreement with that of 
Kimura et al., 14 obtained from macroscopic measure- 
ments. The authors measured a higher critical temper- 
ature at 9 T than ours and they still observed a transi- 
tion at 9.5 T, at about 450 mK, while we did not. This 
discrepancy may arise from a small misalignment of the 
magnetic field with respect to the chain axis in their ex- 
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pcrimcnt, the effect of a magnetic field perpendicular to 
the chain being much weaker. The zoom around the Neel- 
LSDW phase transition (shown in the inset of Fig. 11) 
also evidences the hysteresis of the transition, as already 
mentioned, and, surprisingly, a strong curvature of the 
transition line toward higher field at very low temper- 
atures (below 400 mK). Also, note that the amplitude 
of the hysteresis is getting larger as the temperature de- 
creases (only 0.05 T around 1 K, about 0.15 T at 50 mK). 
Last, between H c ~ 3.9 T and 8.75 T, the agreement 
between our data and the theoretical prediction given 
in Ref. 15 (see dashed blue line in Fig. 11) is excellent, 
while our experimental results do not follow the theory 
any more at higher field. Namely, the system leaves the 
LSDW phase above H p ~ 9.25 T. The NMR measure- 
ments performed by Klanjsek et al.~' confirm our results 
with the same critical field H p and the same transition 
line for the LSDW phase. 



E. Magnetic structure in the LSDW phase 
(H = 4.2 T) 

The magnetic field was set to 4.2 T, coming down 
from 6 T, at T = 50 mK, in order to characterize the 
LSDW phase. These conditions were chosen since they 
correspond to the best compromise in the LSDW phase 
between the strongest possible magnetic intensities and 
the satellites being far apart enough so that they do not 
overlap at all along Ql- The magnetic peak at Ql = 
(from the Neel phase) has also a zero intensity at this 
magnetic field value. In order to determine precisely the 
incommensurate modulation £ of the propagation vec- 
tor k^sDW = (1, 0, £) under these conditions, Ql— scans 
were performed across 12 different pairs of satellites (at 
±£), well distributed in the reciprocal space. They were 
then all fitted to extract the position along c* of the peak 
and their average gave £ = 0.089(4). The magnetic satel- 
lites of the LSDW phase thus appear, under these condi- 
tions, at Q = H ± k LSD w = (h ± 1, k, l± 0.089), with 
h, k, I integers satisfying h + k + I = In. 

The nuclear structure was first checked to be identical 
to that in zero field by performing the nuclear structure 
refinement on a collection of 266 allowed nuclear reflec- 
tions, reducing to 70 independent ones. This refinement 
(agreement factor Rp2 w = 7.72%) also provided us with 
the necessary information for the magnetic structure re- 
finement (i.e., the scale factor, the Co coordinates and 
Debye- Waller factor, the extinction parameters, and the 
A/2 ratio). 34 

For the magnetic structure, 129 different magnetic re- 
flections, reducing to 49 independent ones, were col- 
lected. The magnetic structure refinement was performed 
by continuity to what was obtained in zero field, that 
is: (i) the same structure than in zero field (antiferro- 
magnctic chains with the magnetic moments parallel to 
the chain c— axis and the arrangement between chains 
as shown in Fig. 4) but with an additional modulation 
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FIG. 12: (Color online) Graphical representation of the mag- 
netic structure refinement in the LSDW phase of BaCo2 V2O8 
at T — 1.8 K and H = 4.2 T in decreasing field (coming from 
6 T): calculated versus measured integrated intensities. 



(sine wave function) of the moment amplitude along the 
chains, (ii) the presence of two domains, corresponding 
to each other by the 4i screw axis, as described for the 
Neel magnetic structure. In a first refinement, two pa- 
rameters were fitted, the amplitude A of the sine wave 
function and the domain population, and the following 
was found: A = 1.395(6)/^ with populations 38.6(8)% 
and 61.4(8)%, with the following agreement R— factors: 
R F = 7.97% and R F 2 W = 15.6%. 34 ' 35 Note that the max- 
imum magnetic moment in the LSDW phase is smaller 
than the ~ 2.3/is of the Neel phase. The plot of the 
calculated intensities as a function of the observed ones 
(see Fig. 12) emphasizes the quality of the refinement. 

Figure 13 shows the magnetic structure for domain 
jf=l (using the same labeling as for the zero- field mag- 
netic structure), projected on the (b, c) plane over about 
one half period along the chain axis (~ 5c). The non- 
cquirepartition of the domain populations arising from 
the refinement seems surprising and is discussed in sec- 
tion IV. A second refinement was then performed, im- 
posing equal domain populations (50% each) and the 
following was found: a very similar moment amplitude 
A = 1.364(8)/ig with (as expected) worse R— factors: 
Rp = 13.4% and Rp2 w = 25.8%. As a result, the correct 
analysis of the domain population is not crucial for the 
magnetic moment amplitude, since very similar moment 
amplitude values are obtained in the two refinements. In 
conclusion, this magnetic structure refinement confirms 
for the first time that: (i) the magnetic moment ampli- 
tudes of the four cobalt atoms along the chain (at z =1/8, 
3/8, 5/8, and 7/8) follow the global sine wave function 
imposed by the incommensurate modulation £ = 0.089 
of the propagation vector, thus the magnetic ordering is 
indeed a spin density wave, (ii) the magnetic moments 
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FIG. 13: (color online) Magnetic structure of BaCo2V208 in the LSDW phase at H = 4.2 T (decreased from 6 T) for the 
magnetic domain #1, shown over about one half period (~ 5c), using the same labeling and the same colors as in Fig. 4 for the 
two types of cobalt chains. This projection on the (b, c) plane is to be compared to Fig. 4(a), sketching the magnetic structure 
of domain #1 in zero field. The magnetic domain #2 of the LSDW phase can simply be obtained from domain #1 by reverting 
all spins plotted in red (like for the structure in zero magnetic field). 



are lying along the c— direction, hence this spin density 
wave is longitudinal. 
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FIG. 14: (Color online) Count on top of the 200 reflection 
as a function of the magnetic field evidencing the ferromag- 
netic magnetization M z arising when the field increases. The 
dashed horizontal black line corresponds to the intensity To on 
top of the pure nuclear reflection. The solid black curve corre- 
sponds to the square of the magnetization curve M z of Ref. 1 -1 
corrected for a residual Van Vleck contribution, !i scaled to 
our data, and with an offset To- Inset: ui— scan on the 2 00 
reflection at H = (pure nuclear reflection) and H = 12 T 
(nuclear and ferromagnetic reflection) allowing a precise de- 
termination of the ferromagnetic moment at 12 T. 

In addition to the longitudinal incommensurate com- 
ponent discussed above, a ferromagnetic component is 
induced by the applied magnetic field (field-induced uni- 
form magnetization M 2 ), corresponding to the one ob- 
tained in macroscopic magnetization measurements. The 
latter can also be determined from neutron diffraction, 
as described in the following. The field dependence of 
the square of the uniform magnetization was obtained on 
the complete accessible magnetic field range (0 - 12 T) 



by counting on top of the 200 weak nuclear reflec- 
tion, for which the ferromagnetic component is strong 
(see Fig. 14). Qualitatively, these measurements agree 
perfectly well with the magnetization measurements of 
Kimura et at, 14 that are also represented in the figure af- 
ter subtraction of the residual Van Vleck contribution " 
and after the appropriate treatment for a scaling to our 
data. For a quantitative comparison, rocking curves with 
a long counting time were then performed at 60 mK on 
the 2 reflection at 77 = and H = 12 T (see the inset 
of Fig. 14). Their difference corresponds to the ferromag- 
netic contribution estimated to 0.68/is/Co 2+ using the 
scale factor obtained from structural refinement. This re- 
sult is consistent with the magnetization measurements 
of Kimura et al.: u at T = 1.3 K and H = 12 T, the 
authors measured M z ~ 0.78/iB/Co 2+ yielding M™ 1 ' ~ 
0.69^s/Co 2+ after correction for the residual Van Vleck 
contribution. 33 



IV. DISCUSSION 

The experimental data described in section III reveal 
some features that deserve to be further discussed. Let 
us first discuss the consequences of the existence of mag- 
netic domains, that was probed in this single-crystal neu- 
tron diffraction. The magnetic domains at H = were 
found to be equally populated, as expected. However, a 
sizable non-equircpartition of the domains was found in 
the LSDW phase (populations 38% / 62%) at 4.2 T, as 
mcntionncd in subsection III E. This could be due to the 
shape of the crystal (very elongated along one direction) 
combined with the strong absorption of cobalt. But then, 
the refinement of the zero-field magnetic structure, also 
performed on this crystal, would have yielded the same 
domain populations, which is not the case. 

The non-cquircpartition of the two magnetic domains 
in this field-induced LSDW phase at H = 4.2 T is rather 
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due to a small misalignment of the applied magnetic field 
of the 6 T vertical cryomagnet with respect to the c— axis, 
yielding an asymmctrization of both domains. '' 1 Namely, 
the a— axis was tilted by 1.5° from the horizontal plane 
while the b— axis was tilted only by 0.6°. As a conse- 
quence, the applied vertical field was at 88.5° degrees 
from the a— axis and at 89.4° from the 6— axis, instead 
of 90° from both of them, yielding a nearly three times 
larger field component along the a— axis than along the 
6-axis (0.11 T / 0.04 T, for a 4.2 T applied vertical field), 
and thus favoring one domain with respect to the other. 
Such a large sensitivity of the domain populations to the 
field misalignment was also observed in the magnetoelec- 
tric compound M11PS3. 37 

This non-equirepartition of the magnetic domains 
upon a sizable magnetic field could also explain the be- 
haviors of the field scans in Fig. 7 (performed in the same 
experiment) and Fig. 10 (performed with the 12 T verti- 
cal cryomagnet), as developed hereafter. The field depen- 
dence of the intensity of the magnetic 2 3 reflection sug- 
gests that the populations of both magnetic domains do 
not vary at all when the sweep rate is high [sec Fig. 7(a)], 
while they do vary in the ~ 1 T wide range below the 
transition for a low sweep rate [see Fig. 7(b)]. Indeed, let 
us recall that this reflection is sensitive only to the first 
magnetic domain (see Table III). The only other possi- 
ble explanation for the field scans of Fig. 7(b) would be 
a sudden variation of the magnetic moment in the 1 T 
range below the transition, which seems unrealistic and 
which was never predicted. The small bump observed in 
this field range for decreasing H is also present in Fig. 10 
that was obtained on the same magnetic reflection at an 
even slower sweep rate. Indeed, the tilt of the sample, 
in this experiment using the 12 T cryomagnet, was very 
similar (1.2° for the a— axis, 0.5° for the b— axis). 

Next, we discuss the experimentally obtained critical 
exponents of the order parameter to (i.e., the antifer- 
romagnetic moment) for the transition from the Neel 
phase to the paramagnetic phase above 1.85 K. As this 
is a second-order phase transition, the critical exponents 
characterizing the transitions in the temperature direc- 
tion and in the field direction are expected to be the same 
away from the multicritical points. " s For our case, this 
can easily be shown if we write 



i(H,T) 



to 



1 - 



T 



T C (H) 



(0) 



where T C (H) is the shape of the phase boundary, ft is the 
critical exponent in the temperature direction and too is 
the zero-temperature order parameter. We assume that 
both ft and too are field- independent, which is compatible 
with the experimental results presented in Fig. 5. Taking 
into account that the T C (H) dependence can be locally, 
for small variations of H, approximated by a linear func- 
tion, it is easy to show analytically that the critical expo- 
nent ft' determined by varying H at constant T must be 
the same, ft' = ft. In contrast to this, we experimentally 
found ft' to amount to the value of ft only in the upper 



part of the boundary, at 5 K, whereas we found it to 
decrease in a roughly linear fashion with decreasing tem- 
perature, extrapolating to zero at zero temperature, as 
summarized in the inset of Fig. 6. Such a behavior may 
be an artifact of the method we used to experimentally 
extract the critical exponent ft'. Namely, the method 
that combines fitting to a power law in a progressively 
shrinking field range with the extrapolation to the phase 
boundary, leads to the correct value of the critical expo- 
nent only when ft' depends on the width of the selected 
field range in a linear way. If the exponent ft' reaches 
the value of ft only very close to the phase boundary, in 
the range where experimental fits are no longer feasible, 
the method apparently yields a wrong value. As shown 
in the Appendix, this is exactly the case, so that the ob- 
served temperature dependence of the extracted ft' is a 
mere experimental artifact. 

Now, we wish to discuss the H-T phase diagram we 
obtained and to compare it with the theoretical pre- 
dictions, in the light of the recent NMR results from 
Klanjsek et al. First, our refined magnetic structure at 
4.2 T fully confirms for the first time the realization in 
BaCo2V208 of the predicted LSDW phase, in contrast to 
the previous neutron diffraction study of Kimura et al., 
which demonstrated only that the magnetic order is in- 
commensurate. Second, the first-order character of the 
Neel-LSDW phase transition, first proposed by Kimura 
et al. 14 from field hysteresis in their magnetocaloric mea- 
surements, was confirmed by our neutron diffraction mea- 
surements: (i) a sizable field hysteresis of the transition, 

(ii) a sweeping field rate dependence of both the mag- 
netic intensities and the amplitude of the field hysteresis, 

(iii) the coexistence of the Neel and the LSDW orderings 
over a sizable field range. The first-order character of 
the Neel-LSDW phase transition has not been discussed 
theoretically to our knowledge. 

Last, the LSDW-paramagnctic transition line in our 
H-T phase diagram perfectly coincides, up to 8.75 T, 
with that obtained experimentally by Kimura et al. and 
theoretically by Okunishi and Suzuki. Some deviation 
from these works is however observed at higher fields, 
in agreement with the NMR measurements of Klanjsek 
et al. 2 ' From our neutron measurements and this NMR 
study, the critical temperature of the LSDW phase sud- 
denly drops down to zero when the magnetic field is in- 
creased from 8.75 to H p ~ 9.25 T. This drop corresponds 
to the onset of a new type of magnetic ordering, "' not 
predicted by Okunishi and Suzuki, 15 that occurs before 
the expected transverse staggered ordering 1 1 at higher 
fields (H > H* ~ 15 T). The NMR observations in this 
new phase were explained involving a new type of order 
(named ferromagnetic SDW in their work), where the 
chains with LSDW are ferromagnetically coupled both 
along the a— axis and along the b— axis, with a long- 
wavelength amplitude modulation in the (a, b) plane."' 
Describing the interchain exchange couplings, the au- 
thors were able to develop a model that reproduces both 
the LSDW (named columnar SDW in their work) and 
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ferromagnetic SDW orders. Following these NMR obser- 
vations and their interpretation, the ferromagnetic SDW 
phase was searched in our last neutron diffraction exper- 
iment at H = 12 T and T = 50 mK, by scanning the 
reciprocal lattice around positions compatible with this 
phase. Unfortunately, no signal was observed and addi- 
tional measurements arc thus needed. 



V. CONCLUSION 

The field-induced magnetic properties of the Ising-likc 
spin chain system BaCo2V20s have been studied in de- 
tail by means of single-crystal neutron diffraction. Our 
study was performed with a careful analysis of the mag- 
netic domains distribution, allowing a full determina- 
tion of the magnetic arrangements both in the Neel and 
LSDW phases, at H = and 4.2 T, respectively. The 
determination of the magnetic structure in the incom- 
mensurate phase at 4.2 T is the first direct evidence of 
this predicted LSDW ordering in BaCo2V2 0s. The field- 
induced transitions from the Neel to the paramagnetic 
phase and to the LSDW phase at lower temperature were 
studied by following the positions and intensities of cho- 
sen magnetic reflections up to 12 T and down to 50 mK. 
The phase boundary was found to be first-order for the 
Neel-LSDW transition and second-order for the Neel- 
paramagnetic transition, where critical exponents could 
be extracted from the field and temperature scans, con- 
firming that BaCo2V2 08 belongs to the 3D Ising univer- 
sality class. A complete investigation of the H-T phase 
diagram as well as the field-dependence of the magnetic 
propagation vector was performed and found to compare 
perfectly to the TLL predictions up to 8.75 T. A depar- 
ture from the prediction is observed above 8.75 T where 
the LSDW magnetic satellites disappear. This is well be- 
low the expected transition to the transverse staggered 
long-range order, a discrepancy already pointed out by 
NMR studies and calling for further investigations. 
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VII. APPENDIX 

Assuming that Eq. (6) is valid for temperature sweeps 
across the phase boundary, with tuq and /3 field- 
independent, we show that the critical exponent f3' ex- 
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FIG. 15: (Color online) (a) Zoom on the Neel phase of the 
H-T phase diagram. The same symbols as in Fig. 11 are used 
for the experimental data. Phase boundary is very well ap- 
proximated by the shape of the superellipse as described in 
the text, (b) Dependence of the "critical" power-law exponent 
obtained in a fitting procedure as a function of the width of 
the fitting interval, simulated for field sweeps at various con- 
stant temperatures. As the values close to zero are experimen- 
tally inaccessible, a linear extrapolation of the available points 
(open symbols) to zero (red lines) clearly underestimates the 
critical value, with progressively stronger underestimate as 
the temperature is lowered, (c) Temperature dependence of 
the obtained "critical" exponent j3', to be compared with the 
inset of Fig. 6. 



pcrimentally extracted from the field sweeps across the 
boundary is apparently temperature-dependent. In or- 
der to do that, we convert the magnetization given by 
Eq. (6) into a form appropriate for the field sweeps. For 
simplicity we take the phase boundary with the shape 
of the superellipse, T C {H) = T \/l - {T/T ) n , which 
approximates the experimental boundary very well for 
n = 2.32, T = 5.39 K and H = 3.86 T [see Fig. 15(a)]. 
In addition, we switch to the variable measuring the de- 
parture from the boundary in the field direction, h = 
H C (T) - H = H VI - (T/T ) n ~ H, where H C (T) de- 
fines the transition field at a given temperature T. This 
leads to 



m{h 1 T) = niQ 



1 - 



T/T 



1 - 



^1 - (T/T )» - h/H 







(7) 

For the relevant (3 = 0.32 value and a given constant 
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value of T, this /i-depcndcnce of m is numerically fit, by a 
least squares fit, to a "critical" power law, A(h/H C ) B , in 
the chosen interval < h < hm- We plot thus obtained 
critical exponent B as a function of the fitting range in 
dimensionless units, hat/ H C (T), as shown in Fig. 15(b). 
This procedure is identical to the one used in the exper- 
imental determination of the critical exponent j3', where 
/?' is found by a linear extrapolation of B[hat/T C (H)] to 
zero. The main difference is that in this simulation of the 
experimental procedure, due to the absence of any noise 
and distortions of the data, we have a proper access to the 
required hg± — > limit. Indeed, in Fig. 15(b) we see that 
mathematically correct B(0) = /3 = 0.32 limit is reached 
at any temperature, but in a progressively narrower fit 



range h^ t / H C (T) as the temperature is lowered. 

In contrast to that, very narrow h^t/ H C (T) ranges are 
experimentally inaccessible because the fit becomes un- 
stable. Taking into account that in a typical experimen- 
tal fitting we follow B[hfa/ H C (T)] dependence in a range 
from hat/ H C (T) = 1 down to only 0.3, corresponding 
to open symbols in Fig. 15(b), we see that a linear ex- 
trapolation of these points to zero, B exp (0), gives the 
correct value -B exp (0) = j3 only at high temperature. As 
the temperature is lowered, -B e xp(0) is suppressed in a 
monotonous manner, tending to zero at zero tempera- 
ture [see Fig. 15(c)], reproducing closely the experimen- 
tally observed behavior displayed in the inset of Fig. 6. 



* Corresponding author. Electronic address: grenier@ill.fr 

1 T. Giamarchi, C. Riiegg, and O. Tchernyshyov, Nat. Phys. 
4, 198 (2008) and references therein. 

2 H. Kikuchi, Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. 
Tonegawa, K. Okamoto, T. Sakai, T. Kuwai, and H. Ohta, 
Phys. Rev. Lett. 94, 227201 (2005). 

3 F. Aimo, S. Kramer, M. Klanjsek, M. Horvatic, C. 
Berthier, and H. Kikuchi, Phys. Rev. Lett. 102, 127205 
(2009) and references therein. 

4 B. Grenier, J. -P. Boucher, J.-Y. Henry, L.-P. Regnault, and 
T. Ziman, J. Magn. Magn. Mater. 310, 1269 (2007) and 
references therein. 

5 A. Zheludev, S. M. Shapiro, Z. Honda, K. Katsumata, B. 
Grenier, E. Ressouche, L.-P. Regnault, Y. Chen, P. Vorder- 
wisch, H.-J. Mikeska, and A. K. Kolezhuk, Phys. Rev. B 
69, 054414 (2004) and references therein. 

6 A. Abragam and M. H. L. Pryce, Proc. Roy. Soc. A 206, 
173 (1951). 

7 Z. He, D. Fu, T. Kyomen, T. Taniyama, and M. Itoh, 
Chem. Mater. 17, 2924 (2005). 

8 Z. He, T. Taniyama, and M. Itoh, Appl. Phys. Lett. 88, 
132504 (2006). 

9 Z. He, T. Taniyama, T. Kyomen, and M. Itoh, Phys. Rev. 
B 72, 172403 (2005). 

10 S. Kimura, H. Yashiro, M. Hagiwara, K. Okunishi, K. 
Kindo, Z. He, T. Taniyama, and M. Itoh, J. Phys.: Conf. 
Ser. 51, 9910 (2006). 

11 S. Kimura, H. Yashiro, K. Okunishi, M. Hagiwara, Z. He, 
K. Kindo, T. Taniyama, and M. Itoh, Phys. Rev. Lett. 99, 
087602 (2007). 

12 F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 

13 T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford Univ. Press, Oxford (2004). 

14 S. Kimura, T. Takeuchi, K. Okunishi, M. Hagiwara, Z. 
He, K. Kindo, T. Taniyama, and M. Itoh, Phys. Rev. Lett. 

100, 057202 (2008). 

15 K. Okunishi and T. Suzuki, Phys. Rev. B 76, 224411 
(2007). 

16 S. Kimura, M. Matsuda, T. Masuda, S. Hondo, K. Kaneko, 
N. Metoki, M. Hagiwara, T. Takeuchi, K. Okunishi, Z. He, 
K. Kindo, T. Taniyama, and M. Itoh, Phys. Rev. Lett. 

101, 207201 (2008). 

17 Y. Kawasaki, J. L. Gavilano, L. Keller, J. Schefer, N. B. 
Christensen, A. Amato, T. Ohno, Y. Kishimoto, Z. He, Y. 



Ueda, and M. Itoh, Phys. Rev. B 83, 064421 (2011). 

18 R. Wichmann and Hk. Miiller-Buschbaum, Z. Anorg. Allg. 
Chem. 532, 153 (1986). 

19 P. Lejay, E. Canevet, S. K. Srivastava, B. Grenier, M. 
Klanjsek, and C. Berthier, J. Cryst. Growth 317, 128 
(2011). 

20 J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 

21 Note that for the vanadium atom, the x coordinate (not 
fixed on the 16e Wyckoff site) and the isotropic Debye- 
Waller factor were not refined, since vanadium is purely 
incoherent for neutron. 

22 The phenomenological model for single-crystal refinements 
was used, involving 6 extinction parameters along the crys- 
tallographic axes and their bisectors. But, owing to the 
tetragonal symmetry of the space group and to the crys- 
tal shape (same dimension along a and 6), the follow- 
ing constraints were used for the extinction parameters: 
Ext a = Extb and Ext a +c = Extb+ C - 

23 P. J. Brown and J. C. Matthewman, The Cambridge Crys- 
tallographic Subroutine Library, RL-81-063 (1981). 

24 J. Rossat-Mignod, in: K. Skold and D. L. Price, Editors, 
Methods of Experimental Physics, vol. 23, Academic Press, 
San Diego, CA (1987), pp. 69-157, chapter 20. 

25 Kawasaki et al. 1 ' labeled the propagation vector k,4F = 
(0, 0, 1) which is not compatible with the relation a > c, 
according to Ref. 24, but this does not affect at all the 
result of the magnetic refinement. 

26 E. F. Bertaut, Acta Cryst. A 24, 217 (1968) ; J. Phys. 
Colloques 32 Cl-462 (1971) ; J. Magn. Magn. Mater. 24, 
267 (1981). 

27 M. Klanjsek, M. Horvatic, C. Berthier, H. Mayaffre, E. 
Canevet, B. Grenier, P. Lejay, and E. Orignac, arXiv: 
1202.6376 (2012). 

28 Note that the existence of two satellites in the LSDW 
phase, one on each side of the commensurate position of 
the antiferromagnetic phase, is not related to the existence 
of magnetic domains, as assumed by Kimura et al.: it is 
typical of propagation vectors different from and not on 
a zone boundary of the Brillouin zone. 

29 The extinction was checked to have a negligible variation 
across the transition, so that it does not affect the temper- 
ature and field scans of the magnetic reflections presented 
in Figs. 5 and 6. 

30 V. O. Garlea, A. Zheludev, K. Habicht, M. Meissner, B. 



17 



Grenier, L. P. Regnault, E. Ressouche, Phys. Rev. B 79, 
060404(R) (2009). 

31 A. Pelissetto and E. Vicari, Physics Reports 368, 549 
(2002) and references therein. 

32 M. Mansson, K. Prsa, J. Sugiyama, H. Nozaki, A. Amato, 
K. Omura, S. Kimura, and M. Hagiwara, Physics Procedia 
30, 146 (2012). 

33 The magnetization presented in Fig. 2 of Ref. 14, though 
already corrected for the paramagnetic Van Vleck contri- 
bution by the authors, still exhibits a linear contribution 
XwH of comparable slope below H c and above H a . A fit 
to these data below H c yields a residual Van Vleck contri- 
bution xvv = 0.00733 hb/T. 

34 Note that the agreement is not as good as for the cubic 
crystal in zero field (see subsections IIIB for the nuclear 
structure, IIIC for the magnetic structure), which is due 



to the shape of this crystal (elliptical instead of cubic), for 
which the absorption is more difficult to correct. 
An additional reason for the poorer Rf— factor of the 
LSDW phase refinement as compared to the Neel phase 
are the 7 times smaller intensities of the magnetic reflec- 
tions in the former. 

Using the UB matrix obtained from the centering of about 
30 reflections in the experiment using the 6 T vertical cry- 
omagnet (in which the LSDW magnetic structure was re- 
fined and the field scans of Fig. 7 were measured), the 
precise orientation of the a and b axes with respect to the 
vertical magnetic field could be calculated. 
E. Ressouche, M. Loire, V. Simonet, R. Ballou, A. Stu- 
nault, and A. Wildes, Phys. Rev. B 82, 100408(R) (2010). 
W. Paul, J. Magn. Magn. Mat. 87, 38 (1990). 



